Dirac fields in curved spacetime as Fermi-Hubbard model with non unitary tunnelings 
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In this letter we show that a Dirac Hamiltonian in a curved background spacetime can be in- 
terpreted, when discretized, as a tight binding Fermi-Hubbard model with non unitary tunnelings. 
We find the form of the nonunitary tunneling matrices in terms of the metric tensor. In a simple 
case of a static diagonal metric, the tunnelings become unitary. Alternative interpretation of the 
Fermi-Hubbard Hamiltonian is that of a Pauli Hamiltonian, i.e. a non relativistic limit of the Dirac 
Hamiltonian. In this case the tunnelings remain, in general, non unitary even for the static diag- 
onal metric. We discuss a possibility of synthesizing such Hamiltonians by means of laser assisted 
tunnelings in cold atomic experiments. 
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Recently, much interest was devoted to a study of many 
body physics of quantum gases [l|, |2j. High degree of 
control of the experimental parameters has allowed for 
designing specific Hamiltonians [3(- A special category 
is then a fabrication of synthetic gauge fields [4| , where 
a remarkable experimental progress has been achieved in 
last couple of years, including the realization of synthetic 
electric [fj and magnetic [6[ fields in the bulk as well as 
on the lattice [7| . A non abelian synthetic gauge field of 
the spin orbit Rashba type has been demonstrated in the 
bulk |8J, lS| . In the case of a lattice, which will be of main 
interest, laser assisted hoppings lfj,|ll| allow for a sim- 
ulation of a (non abelian) lattice gauge theory [l2| with 
cold atoms [l3|. Different works adressed the question 
of non abelian background fields with cold atoms. Such 
situation occurs e.g. in the case of electrons with spin or- 
bit coupling, in both non interacting |14| and interacting 
[15| cases. In those scenarios, however, the tunnelings of 
different spin components between two adjacent lattice 
sites are described by unitary matrices, due to the her- 
miticity of the gauge field [12|, llfj. Moreover, an explicit 
form of these matrices is determined by the theory one 
wants to simulate, e.g. the mentioned spin orbit cou- 
pling. It is thus an interesting question, what happens, 
if these tunnelings become non unitary, which can be, 
in principle, done in the cold atomic experiments using 
known techniques, as explained below. 

Driven by this motivation, let us start with a kinetic 
fcrmionic Hamiltonian in d — 1 spatial dimensions of a 
form 



nates as x (= (t, x)), while the space coordinates as x. 
The matrices T(x, x + ak) represent a parallel transporter 
of a quantum field ip between sites x and x + ak ■ In the 
case of lattice gauge theories, the matrix T belongs to a 
representation of a gauge group, which is typically com- 
pact, such as U(n) with n being the number of "flavor" 
components of the field vb. In such case, the matrices T 
become unitary. The relevant question is thus, what if 
the T are non unitary? This question has actually been 
already adressed in the past [17H19J and more recently 
[2 fj| | , but such interpretation seems to be problematic and 
was not actively pursued. Lets return to the Hamiltonian 
Eq.([T]). In what follows, we will be interested in physics 
of fermions, so that ip is a spinor. As discussed later, a 
general matrix T <= Gl(n,C) can be engineered in cold 
atomic systems (n is the number of spin components). 
We would like to emphasize, that for such a novel situ- 
ation, the Hamiltonian Eq.JTJ) is interesting on its own 
right. However, it is interesting to look, whether some 
physical significance can be given to it. 

A starting point of our discussion will be a classical 
(field not quantized) fcrmionic field in a curved space- 
time, which can be described by a Lagrangian density 

El 



C(x) = *Jg \ -ivj{x)Y i D^vj{x) +h.c. - mip{x)vb(x) 

'(3) 

Let us recall [22j, that working in a coordinate basis e M , 
in which the spacetime vector x is defined in terms of 



H{t)= Y ^l( x )T s ^(x,x+a k )^ s ,( x +a k )+h.c, (1) its components a?/*, x = x%, one may construct a local 



where the sum runs over all lattice sites x and directions 
k = l..d— 1 and the fermionic operators satisfy the usual 
anticommutation relations 



{* s (x),^,(x')}=W<W. (2) 

Throughout the paper, we denote the spacetime coordi- 



orthonormal basis e Q . The two bases arc related through 
vielbeins e a = e^e^. The metric tensor g is defined as 
Va/3 — Ra e 09fii', where r\ is the Minkowski, i.e. flat metric. 
Then, the curved spacetime 7 matrices are defined as [23[ 
7 M = 7 a e^, where 7" are the usual (flat spacetime) Dirac 
matrices, for which {7°, 7 /3 } = r\ a ^ and we adopt the sign 
convention 77 = (+,—,—,...). The covariant derivative 
acting on the spinor is D^ = d^ — T^, where r^(a;) = 



i [ 7 a , 7 ' 3 ] e£(x)(V ^epv) and V M e^ = d^ep u - T^e^ 
(for a brief overview with essential technical details, see 
[24|). The canonically conjugate momentum to if) can be 
found in a usual way as 



" (x) = afM = ^2*° 



(4) 



and similarly for 7r which is conjugate to if). One then 
obtains the Hamiltonian density H = ir(doip) + (doif))ir — 
C, 

H{x) = -^gl^iip [x k D k - 7°r ] ip + h.c. - mipip 

'(5) 

Lets now cosider an isotropic square lattice in coordinate 
basis with lattice spacing a. We introduce a covariant 
derivative on the lattice as 

D^(x) = - [P(x, x + a^ipix + a M ) - 4>(x)} , (6) 

where P{x, x + a k ) is the parallel transporter from x + a k 
to x and reads (here fi is fixed) 



P(x,x 



"Pcxp 



darT^ar) 



X + CLu 



(7) 



and V stands for the path ordering. One can then for- 
mally discretize the Hamiltonian H(t) ~ J d d ~ 1 xH(x), 

Hit) = ^Yj ^(x)M k P(x, x + a k )if>(x + o fc ) 



x.A- 



if)\x) 



aM T + -^gm-f 



ip(x)+h.c, (8) 



where M^{x) = — i-Jg^f ^. At first sight, the structure 
of the Hamiltonian Eq.© is similar to Eq.©, but there 
are two major differences. First, in the former case, the 
fields are not quantized and second, they are time de- 
pendent, so it is not obvious, whether one can obtain the 
desired anticommutation relation Eq.© for the space de- 
pendent operators. In order to proceed, we shall rely on 
the arguments exposed in [25j. We summarize the main 
steps crucial for our purpose. The classical field ip{x) is 
a projection of ket \ip) to a spatiotcmporal basis \x) 



if)(x) = (x\ip) 



(9) 



Similarly, one obtains a different field which is only space 
dependent on some constant time hypersurface as 



vfr(x) 



(10) 



The relationship between the two fields if)(x) and \&(x) 
can be found from the equivalence {(f>\ip) = (4>, ip), where 
the scalar product in curved spacetime is defined as [25J 



,V)= / d'-V^VYV- 



(ii) 



From the resolution of identity J d d 1 cc \x) v / 7?7°7° (x\ — 
1 one obtains the evolution of the ket \x) 

do \x) = -\ \x) (doVgj l )(V9J°l )-\ (12) 
which can be formally integrated to give 

|x) = |x,t) = V2|x)(V7i7 7 )^- (13) 

The meaning of the factor y/2 will become clear momen- 
tarily. 

The quantization of the space and time dependent field 
Eq.©, ip(x), then proceeds by imposing equal time anti- 
commutation relations for the canonically conjugate op- 
erators [23[, namely 



{if> s {x),TT s ,(x')} = iS(x- x')S s 



(14) 



where ir(x) is given by Eq.©. We can use the relations 
Eq. © [lOl) and Eq. (fT3|) to find the relationship between 
the two fields ip(x) and \I/(x) to be 



V(a;)-V2(V57 7 ) 5 *(x). 



(15) 



Plugging Eq.© into Ea. ([T4")) and using Ea. (Ti"5|) we obtain 
for the anticommutator of the constant time hypersurface 
fields 

{vMx),^,(x')}=<5(x-x')<W, 

which is precisely the relation Eq.© (the factor y/2 in 
Eq. ([13[) is actually an integration constant determined 
by the form of the anticommutator Eq.© and originates 
in the hermitian definition of the Lagrangian Eq. © con- 
taining factors 1/2). 

Lets now take the lattice Dirac Hamiltonian Eq.© and 
write it as 

H(t) = J2^(x)T(x 7 x+a k )ip(x+a k )+h.c.+^(x)V(x)ip(x), 

X./v 

(16) 
We now use the relation Ea. (fT5"]) to substitute for if){x) 
and write the Hamiltonian as 

H(t) = Y,*Hx)T(x,x+a k )y(x+a k )+h.c.+^(x.)V(x)y(x.), 

X,fe 

(17) 
where 

T(x,x + a,i) = 

^(iM°)- 1 / 2 {x)M l P(x, x + a^diMyy 1 / 2 ^ + a,) 

V(x) = 

(zAf°)- 1/2 [ M °r + (M°r )t + Jgmj ] ((*M°)t)-^) 

where we have put the lattice spacing a = 1 for simplic- 
ity. Now, the Hamiltonian Eq. (fTT)) has the same struc- 
ture with the correct anticommutation relations for the 



operators as Eq.flTJ (plus the local term). The price to 
pay in order to achieve this goal was to absorb the spa- 
tiotemporal dependence of the fields tp(x) to the elements 
of the Hamiltonian and thus spoiling its covariance. 

Next, we discuss a non relativistic limit of the Dirac 
equation. After all, the kinetic part of the usual Hub- 
bard model for electrons in tight binding approximation 
(Eq.(TT]) with T ex 1) is obtained from the non relativis- 
tic quantum mechanical Hamiltonian H ex p 2 /(2m). It 
is thus interesting to see, how similar derivation works 
for fermions in curved spacctimc background. We should 
use a systematic method, known as Foldy-Wouthouyscn 
transformation [26j (also used in the context of quantum 
fields in curved spacetimes J27| ) , which perturbatively de- 
couples the electron and positron modes. One can derive 
the Dirac equation from Eq.Q 

(iYD^ - m)i){x) = 0, (19) 

which can be rewritten as Schrodingcr equation 

td ^x) = n D ^(x), (20) 

where 

H D = ( 1 °)- 1 (rn-i 1 k D k ) + ir (21) 

is the Dirac Hamiltonian. The non relativistic limit can 
be obtained from the Dirac Hamiltonian of the form 



H D = l°m + £ + 0, 



(22) 



where £ and O are even and odd operator, defined by the 
property [7°,f] = and {7 ,O} = and 7 is in the 
Dirac representation. The lowest order expression for the 
nonrelativistic Hamiltonian is 



1 



2m 



(23) 



where the subscript P stands for the Pauli Hamiltonian. 
One can identify 7 , £ and O by comparing Ea. (j2"2|) with 
Eq. ([2"!]) . In the most general case it yields rather lengthy 
expressions. In order to proceed with the calculation, we 
will thus consider a simple, yet non-trivial scenario with 
a static diagonal metric of the form 



1 
h 



and h = dia,g(ha(x k )), where i = l..d— 1 and the diagonal 
terms depend only on spatial coordinates x k . First thing 
we note, is that in this case, the viclbcin fields are also 
diagonal, e^ = for /i 7^ a. In particular e[J = 1 implying 
(7°) _1 = 7° = 7°- Also > r o = 0- We then obtain for the 
Dirac Hamiltonian 



H 



D 



7 m 



i-f°j k D k = 7°tb + O, 



since the term 7°7 fe .D/ c is odd for the metric considered. 
We then obtain for the Pauli Hamiltonian density 

Up = 7 °m - -L 7 (7 7^ fc )(7°7 J ^)- (26) 
2m — — 

The total Hamiltonian, expressed in terms of field vari- 
ables, then reads 

H P = ^[(ij,Hp^) + (Hpij^)}. (27) 

The scalar product can be evaluated by integrating per 
parts in curved spacetime. The reason why one wants to 
do that is to obtain terms of type (Dip' ! )(Dilj) rather than 
ip^D 2 ip, since the former can be mapped to a Hubbard 
luodel with only nearest neighbor hopping. Evaluating 
Eq.(H7]), we get 



H f 



1 1 

22m 



id-l 



x^/g^ k i° D k D j ip+h.c- 



id-l 



Xy/gmiptp 
(28) 



At this point ip is still 2l d / 2 l component spinor, where [nj 
is the integer part of n. By construction, the Hamiltonian 
Hp contains only even operators and we can thus split 
the spinor into two parts, say tp = (x,tp)- Each of the 
spinors x, <p has 2^ d ' 2 '~ 1 components, which will have in- 
dependent dynamics. In case of the diagonal static metric 
and d = 4, we find 



tViw* = - ( e * e ^ v ^ 



'ib-^j 







( 29 ) 
- Tfe, Tfc = 

Lets write the Pauli 



where Vfe = d k 

Hamiltonian for one of the spinor components, say \, 
which we write as 



/ dV/fiViViX + j d^" 1 xVfcicrV^ [V fe , V,-] \ k<j m 

where f k j = Jlje\&\- The commutator in the second 
term is familiar from non-abclian gauge theories and we 



(24) h a ye[V k ,y j ] = d [j T k] 



r fe ,r, 



which acts locally on 



(25) 



the spinor x- The first term can be integrated per parts 
to yield (using ft = — f,) 

dV/nViViX = 

d 3 x {MV iX )HViX) - (difa)xHViX)} ■ (31) 
We thus write the Pauli Hamiltonian as 



d 3 x- — 
2 2m 



/«(V lX ) t (Va) - ($/«)** (Vtf) - X%^MV fc ,V,]L . X ,l +li.c.+VsmA 



(32) 



r 



We are now in the position to discretize the Pauli Hamil- 
tonian, which is to follow exactly the same steps as in 
the case of Dirac Hamiltonian. Using again the prescrip- 
tion Eq. (|15[) . which now takes a simple form x( x ) — 
\/2y/g~2X(x.), we arrive at a Hamiltonian, which can 
be formally written as Eq. (|17[) , where wc have to replace 
"L — > X and the matrices T, V now depends only on spa- 
tial coordinates x and read 



(>,) 



T(x,x 
V(x) = 2m 
V9- 1 



V9~ 



-(/« 



1 



ft/ il )P(x,x + a i ) 



fu + fu-(difu) + 



jVssVlVfc.Vi] 



fc<J 



+ h.c.( 



where /^~ = /ji(x — a^) and we have to replace T — > T in 
the definition of the parallel propagator Eq. © . 

Physical interpretation; we just provided two possible 
interpretations of the kinetic Hamiltonian of the form 
Eq.([17p. driven by a possible interpretation of non uni- 
tary T. However, it is interesting to notice, that for the 
static diagonal metric Eq. (|24p . 7 are unitary, T^ are anti- 
hermitian, Tj, = —T^ and the parallel propagators Eq.([7]) 
become unitary. We thus have T which is also unitary, 
contrary to the Pauli Hamiltonian case. Another ques- 
tion is what field theory we are actually simulating. Lets 
take an example of simulation, where the Hamiltonian 
Eq.([T]) describes a motion in a (two dimensional) plane 
for a two component field xj} s , s = 1,2. If we want to 
interpret it as a Dirac Hamiltonian, our simulator would 
correspond to a Dirac Hamiltonian in 2+1 dimensions. If 
we are to interpret it, however, as a Pauli Hamiltonian, 
the simulator corresponds to a spin half fcrmion living in 
3+1 dimensions, but whose motion is confined to a plane. 

Implementation with cold atoms; in the experiments 
with cold atoms, the internal degrees of freedom are usu- 
ally played by the hyperfine states of the atoms. These al- 
low for a laser assisted tunnelings between adjacent sites, 
say i, i 1 of the optical lattice. Lets denote the internal 
degrees of freedom s. In order to engineer an arbitrary 
T(x) £ GZ(2,C), it is necessary to control each of the 
tunneling rates (i,s) -f-> (i',s') independently in all spa- 
tial directions and moreover, the rates in general vary in 
spacetime. Different techniques and their combination 
can be used in order to achieve this goal. For example, 
bichromatic lattices can be combined with an indepen- 
dent Raman laser for each transition s -H- s' [28[. The 
spatial dependence is then given by a transverse profile of 
each Raman laser. It can be given e.g. by a (typically) 
gaussian laser profile which varies slowly on the lattice 
spacing or it can be designed using ejj. a specific phase 
masks |29f or array of microlenses [3(3 ] . which allow for 



the modulation on the scale of lattice spacing and were 
already used in the cold atomic experiments. Another 
comment is, that the potential V{x) in Eq. (|17[) is non 
diagonal and might be difficult to engineer. The way 
around is that since V is hcrmitian, it can be diagonal- 
ized by unitary transformation. It amounts to redefine 
the tunneling matrix T (analogous to a local gauge trans- 
formation in the case of gauge fields) and the spinors \l/. 
Since the transformation is unitary, the anticommutation 
(3ft)lations Eq.® are preserved. 

Interactions; so far, we were talking about non inter- 
acting case. Although non abelian lattice gauge theories 
Ke non trivial already at this level, the most interesting 
physics can be obtained in the presence of interactions. 
A natural interaction term for spin half fcrmions in op- 



tical lattice is H irA <xUY^ 



S^S' 



n s n s 



where n s 



iptfa 



the density operator. Once again, one entirely legitimate 
approach is to consider a Hamiltonian H = H\d n + -Hint , 
with H kin Eq.© and T £ Gl(2,C) as such and study 
its properties (H could also describe interacting bosons 
instead of fermions) . The other approach is to design it 
in a way, that it simulates a given field theory. For ex- 
ample, a proposal of simulation of a Thirring model (i.e. 
1+1 dimensional field theory) with cold atoms was made 
recently [31[, where the interaction term reads J^J^ with 
J^ = ipj^Tp. In curved spacetime, the replacement 
7 M — > 7^ makes the interaction term spacetime depen- 
dent. One can thus try to modify the proposal |3l| in 
a way that creates the correct interaction term, which 
might be an interesting test bed situation, since as one 
dimensional theory, the massless Thirring model is solu- 
ble also in curved spacetime 32 [. 

We would like to mention, that a simulation of a Dirac 
field in curved spacetime with cold atoms was already 
331 ] . but the discretization was carried out in 



adrcsscd 

the limit of small lattice spacing, such that the approx- 
imation P w 1 + aT is valid. For stationary metrics, 
considered in [33J, it results in unitary tunnelings. 
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